1 



LOCAL BENCHMARKS OF INFRARED GALAXY EVOLUTION : THE SWIRE-SDSS DATABASE, 
FAR-INFRARED LOCAL LUMINOSITY FUNCTIONS & VIRTUAL OBSERVATORY TOOLS 



Mattia Vaccari', Lucia Marchetti''^, Alberto Franceschini', and Ismael Perez-Fournon^ 

^Department of Astronomy, University of Padova, Italy, inattia@inattiavaccari.net 

^Instituto de Astrofisica de Canarias, La Laguna, Tenerife, Spain 



ABSTRACT 



We describe the construction and the properties of the 
SWIRE-SDSS database, a preliminary derivation of the 
Far-Infrared Local Luminosity Functions at 24/70/160 
micron based on such a database and ways in which VO 
tools will allow to refine and extend such work. 
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1. INFRARED GALAXY SURVEYS ACROSS 
COSMIC TIME 

In the context of ever deeper surveys at most wave- 
lengths, it is becoming increasingly difficult and impor- 
tant to reliably measure galaxy properties in the Local 
Universe: difficult because the very possibility to carry 
out extremely deep observations leads to most observ- 
ing time being spent on the deepest pencil-beam surveys 
rather than on shallower wider-area ones which are grad- 
ually becoming the domain of dedicated facilities (see 
e.g. 2MASS, SDSS, Pan-STARRS and LSST) and impor- 
tant because the increasingly detailed knowledge of the 
high-redshift Universe needs similarly well-defined local 
benchmarks to trace the formation and evolution of galax- 
ies across cosmic time in great detail. Perhaps more im- 
portantly, in the era of multi-wavelength surveys and vir- 
tual observatories, shallow wide-area surveys with large 
data rates are likely to profit the most from the paradigm 
shift caused in astronomical research by the easy access 
to a number of otherwise separate databases for science 
exploitation. 

This is particularly relevant for infrared galaxy surveys 
because, as first demonstrated by IRAS, infrared galax- 
ies undergo a dramatic evolution in the Local Universe. 
As later shown by ISO and Spitzer, such an evolution is 
sustained up to z ^ 1, when it subsides leading to an 
approximately constant infrared luminosity density up to 
z ~ 3. Accurately pinpointing the infrared local luminos- 
ity function and studying its evolution over the < z < 1 



range is therefore of the uttermost importance for studies 
of galaxy formation and evolution. 

This work capitalizes on the above trends by exploit- 
ing (mainly) two major survey projects, in the infrared 
and optical respectively, which have recently got close to 
completion. The SWIRE (infrared) and the SDSS (opti- 
cal) datasets in the Lockman Hole region, along with an- 
cillary datasets at optical and near-infrared wavelengths 
such as 2MASS and UKIDSS, are used to derive the 
galaxy local luminosity function at 24, 70 and 160 /im 
and thus place stronger constraints on models for the for- 
mation and evolution of infrared galaxies. 



2. THE SWIRE-SDSS DATABASE 



The SWIRE (Lonsdale et al. 2003) and SDSS (York 
et al. 2002) projects provide remarkably well-matched 
datasets for studies of optical vs infrared galaxy emission. 
This is particularly true in the Local Universe, where 
redshifts were measured spectroscopically by the SDSS 
for the majority of infrared- (as well as optically-) bright 
sources. 

The size of the catalogs produced by these projects 
at multiple wavelengths, however, requires robust auto- 
mated procedures to deal with the assessment of cata- 
log quality (and most importantly its reliability and com- 
pleteness), the association of sources detected at differ- 
ent wavelengths and the optimal handling of multiwave- 
length photometric and spectroscopic datasets, most of 
which coming with their own value-added data products 
such as photometric and spectroscopic redshifts, stellar 
mass and star formation rate estimates etc. 

For our purposes, i.e. in order to determine the LLF 
of FIR-selected, there are two critically important steps 
while handling the data. Namely, one first needs to deter- 
mine the correct optical and NIR counterpart of each FIR 
source and then go on to determine a reliable redshift for 
each source based on either photometric or spectroscopic 
means. 
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The conceptually simple issue of identifying a source 
in images taken at different times, spatial resolutions or 
wavelengths, has recently been receiving a great deal of 
attention, largely due to the increasing source densities 
provided by modern detectors and coordinated surveys 
carried out at a number of wavelengths. The nearest- 
neighbor traditional approach has thus been gradually su- 
perseded by the likelihood ratio technique by Sutherland 
& Saunders 1992 as well as by other more sophisticated 
methods such as the fully Bayesian approach by Budavari 
& Szalay 2008. Similarly, imaging in 5 or more bands 
and/or at NIR wavelengths greatly improve the accuracy 
of photometric redshifts, as demonstrated, between oth- 
ers, by Rowan-Robinson et al. 2008 and Ilbert e al. 2009. 

In our case, limiting the study to bright (>S'24.7o.i60 ^ 
1, 15,75 mJy) and local (0 < z < 0.25) FIR sources 
greatly simplified both issues. On one hand, the source 
density of bright FIR sources is small enough to provide a 
straightforward optical identification, and particularly so 
when suitably deep imaging at intermediate IRAC wave- 
lengths are available. On the other hand, spectroscopic 
redshifts are provided by the SDSS for a large portion of 
bright local galaxies, while carefully calibrated and vali- 
dated photometric redshifts are available for virtually all 
bright sources via the SDSS collaboration. 

In this work we therefore used the upcoming SWIRE 
final data release and the SDSS DR6 in the Lockman 
field, where the joint coverage by SWIRE and SDSS 
amounts to ^ 10 deg^. Nearest-neighbor search be- 
tween MIPS70 & MIPS 160 micron FIR single-band cat- 
alogs and IRACH-MIPS24 MIR 4-band catalogs guaran- 
teed a reliable and complete association process. Simi- 
larly, association between IRACH-MIPS24 MIR catalogs 
and SDSS 5-band catalogs provided a robust link to the 
optical, along with a vast spectroscopic and photometric 
redshift database. A number of ancillary optical (SWIRE 
team proprietary imaging) and NIR catalogs (2MASS 
PSC and UKIDSS DR4) were also positionally corre- 
lated against SDSS, providing additional source diagnos- 
tic power over the full field (albeit with 2MASS PSC 
only) but also deeper photometry over portions of the 
field. Source and redshift reliability flags were extracted 
from the relevant catalogs and included in the database. 
Stars (only a few at the bright FIR flux limits we were 
considering) were identified on the basis of a number of 
criteria, so as to provide an extragalactic source catalog 
of the uttermost quality. 

The basic properties of the SWIRE-SDSS 70 micron 
sample, here selected with a flux cut of 25 mJy, which 
is where SWIRE completeness stars falling off, are illus- 
trated in Figure[T] For such relatively bright FIR sources, 
the optical identification rate (first panel) turns out to be 
very high, with most of the sources which failed to be 
reliably assigned either an optical counterpart or a red- 
shift believed to reside at higher redshifts. Their mag- 
nitude distribution (second panel) sharply falls off well 
before the SDSS magnitude limit, guaranteeing an accu- 
rate magnitude (and thus redshift) measurement for most 
sources. The redshift distribution (third panel), although 



peaking at low redshifts, shows an extended tail, high- 
Ughting the potential of the SWIRE-SDSS FIR-selected 
database for investigating the higher-redshift universe 
and thus the evolution of the luminosity function from 
the local universe up to at least z > 0.5. The fraction of 
spectroscopic redshifts decreases sharply at z ~ 0.15 and 
r ~ 18, as expected on the basis of SDSS spectroscopy 
selection criteria. The fraction of photometric redshifts 
deemed "unreliable" by the SDSS DR6 photometric red- 
shift pipeline is limited throughout most of the redshift 
and magnitude range, increasing sharply only at z > 0.6 
and r > 22. Finally, the luminosity distribution of lo- 
cal sources (fourth panel) spans more than three dex and 
thus allows to place robust constraints on the luminosity 
function. 



3. THE SWIRE-SDSS FIR LLPS 



Armed with the reliable and complete sample of bright 
FIR sources provided by the SWIRE-SDSS database, we 
computed the monochromatic luminosity functions over 
< z < 0.25 at 24, 70 and 160 micron applying the 
^/Vmax estimator to 524,70,160 > 1, 15, 75 mJy sources. 
A completeness correction based on simulations was ap- 
plied to 70 and 160 micron sources fainter than 25 and 
125 mJy respectively. In Figures |2] and [3] we compare 
our results with previous work as well as with a back- 
ward evolution model for the evolution of FIR sources 
by Franceschini et al. (in prep). The black solid lines 
indicate the expected total LF, while the four contribut- 
ing populations are respectively: slowly or non-evolving 
disk galaxies [blue dotted lines]; type-1 Active Galac- 
tic Nuclei evolving as shown by UV and X-ray selected 
quasars and Seyferts [green long-short dashed lines]; 
moderate-luminosity starburst galaxies with peak emis- 
sion at z 1 and type-2 Active Galactic Nuclei [cyan 
dot-dashed lines]; ultra-luminous SCUBA-like starburst 
galaxies with peak evolution between z = 2 and z = 4 
[red long dashed lines]. The agreement with the model is 
mostly very good, albeit within the uncertainties affecting 
the sample particularly at fainter luminosities. Addition- 
ally, our results overcome the large uncertainties charac- 
terizing some previous results. Note that errors in data 
points reported from previous work are not shown for 
clarity but are in most cases larger than ours, and that 
the full SWIRE-SDSS 70 micron sample will eventually 
include more than twice the sources used in this work. 



4. HARNESSING THE POWER OF VO TOOLS 



This work provides further proof of the power of archival 
astronomy in the internet age by combining value-added 
data products by two large-area legacy survey projects 
in the infrared and optical wavelength range. Streamlin- 
ing access and visualization of large datasets as well as 
data processing steps combining data acquired at differ- 
ent wavelength, VO tools have greatly contributed to en- 
abling similar studies of an increased depth and breadth. 
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While the prehminary analysis of the SWIRE-SDSS FIR 
Local Luminosity Functions presented in this work re- 
lied on a set of ad hoc scripts mostly developed in IDL 
and SuperMongo, in order to refine and extend our re- 
sults, we will be applying tools developed within the VO 
to streamline the workflow and make the results available 
to the community in a VO-compliant way. 

An exploratory effort at reproducing the steps we fol- 
lowed in catalog inspection, visualization and cross- 
correlation using TopCat has not only made very rapid 
progress but also enabled much more thorough checks 
on various aspects of data quality thanks to TopCat built- 
in interactive plotting facilities. In the future we plan 
to exploit as much as possible the scripting capabilities 
provided by the VO Desktop and AstroGrid Python, and 
look forward to when these will support IDL scripting as 
well. This will arguably simplify to a large extent the im- 
plementation of a number of improvements, such as the 
extension to SWIRE ENl and EN2 fields (thus almost 
doubling covered area), and the use of SDSS DR7 and 
UKIDSS DR5. 

A reliable and extensible VO platform and the availabil- 
ity of high-quality well-documented datasets via online 
archives will be the main factor driving an increased 
adoption of VO tools by the community in the near fu- 
ture. An outstanding opportunity to get the community 
excited about the VO will be offered by the launch of 
Herschel, an ESA Cornerstone mission observing in the 
FIR & SMM wavelength range and the first "great ob- 
servatory" to be launched in the VO era. Herschel Data 
Processing tools will be tightly integrated with VO tools 
through a Plastic server, offering a seamless way to ex- 
change data in tabular and image formats. This will en- 
sure an optimal early science exploitation of Herschel 
data, and incidentally allow to detail and improve con- 
straints on the FIR LLF estimates presented in this work. 
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Figure 1. The SWIRE-SDSS 70 micron sample. From the 
top : identification rate, magnitude, redshift and luminos- 
itv distribution. 
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Figure 2. SWIRE-SDSS MIPS LLF. 24 micron - Filled Squares (This Work), Empty Squares from Shupe et al. 1998, Empty 
Triagles from Marleau et al. 2007. Lines are models from Franceschini et al. (in prep). See text for details. 




Figure 3. SWIRE-SDSS MIPS LLF. Left Panel - 70 micron - Filled Squares (This Work), Empty Squares from Saunders 
et al. 1990, Empty Triangles from Takeuchi et al. 2003. Right Panel - 160 micron - Filled Squares (This Work), Empty 
Triangles from Takeuchi et al. 2006. Lines are models from Franceschini et al. (in prep). See text for details. 



